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Abstract 



In this paper we show how quantum corrections, although perturbatively 
small, may play an important role in the analysis of the existence of some 
classical models. This, in fact, appears to be the case of static, uniform- 
density models of the interior metric of cosmic strings and neutron stars. We 
consider the fourth order semiclassical equations and first look for perturba- 
tive solutions in the coupling constants a and (5 of the quadratic curvature 
terms in the effective gravitational Lagrangian. We find that there is not a 
consistent solution; neither for strings nor for spherical stars. We then look 
for non-perturbative solutions and find an explicit approximate metric for 
the case of straight cosmic strings. We finally analyse the contribution of the 
non-local terms to the renormalized energy-momentum tensor and the possi- 
bility of this terms to allow for a perturbative solution. We explicitly build up 
a particular renormalized energy-momentum tensor to fulfill that end. These 
state-dependent corrections are found by simple considerations of symmetry, 
conservation law and trace anomaly, and are chosen to compensate for the 
local terms. However, they are not only ad hoc, but have to depend on a 
and (3, what is not expected to first perturbative order. We then conclude 
that non-perturbative solutions are valuable for describing certain physical 
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situations. 
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I. INTRODUCTION 



It is well-known that semiclassical quantization of matter fields in a curved background 
reveals the necessity to include higher order terms in the gravitational Lagrangian in order 
to absorb the divergent behavior coming from the one-loop expansion of the matter field 
(see Ref. [0 for a review). These terms can be expressed as quadratic combinations of the 
curvature tensor and of its derivatives. The resulting theory has field equations including 
fourth order derivatives of the metric. 

The loop expansion to quantum gravity can be formally interpreted as an expansion in 
powers of l 2 { (i. e. in powers of h in G = c = 1 units). One formally expands both a 
quantum metric operator g^ u = r/^ + ^/^ u and a quantum matter field operator <J> Mi , (for 
simplicity a scalar field) about any fixed vacuum classical solution gffj. However, if we 
consider the full series in this expansion we discover that not only does not converge but is 
non-renormalizable. In the lowest (one-loop) order, the matter divergences can be absorbed 
in the constants of the theory. However, this finiteness cannot be attained at this one-loop 
order for gravity coupled to matter fields. This raises the well-known graviton loop problem 
since, there is no obvious mathematical nor physical justification for dropping the 'one-loop' 
of the graviton, because in this scheme the one-loop quantum effects of the metric are just 
as important as those of the matter field. 

On the other hand, the alternative 1/N approach to semiclassical gravity was investigated 
in Ref. |[22|| . In this approach one considers the Einstein-Hilbert gravitational action plus 
matter (|1 . 1|) , where in the matter Lagrangian C 1 ^ one assumes the presence of N identical 
and free (i. e. non-self interacting) conformally invariant scalar fields, all of which are in the 
same quantum state. Then, instead of expanding the theory in powers of h, which results 
in the ordinary loop expansion, one expands in powers of 1/N. In this expansion each field 
is now coupled to gravity with coupling proportional to 1/N and all matter contributions 
(classical and one-loop) are included already at the leading order in 1/N. Note that one of 
the advantages of such an expansion is that it is a manifestly gauge-invariant expansion (N 
being a gauge invariant parameter). It was also observed that it is easier to examine 



the theory in the limit N — > oo, where even the one-loop graviton operator and Faddeev- 
Popov ghosts contributions can be neglected with respect to the contribution coming from 
a large number N of matter fields. Thus, in this approximation there are no corrections 
from the quantized degrees of freedom of the gravitational field itself and gravity can be 
simply treated classically. Moreover, observe also that the N fields, although interacting 
with gravity, do not couple with another possibly present matter field. Note also that in the 
N —>■ oo limit, the fluctuations in the expected stress-energy tensor of the N matter fields 
become negligible. The applicability of the semiclassical field equations, Eqs. (|1.2| ) below, is 
now given by the higher order l/iV-corrections, which suggest that these equations would 
(in principle) be valid even near the Planck scale, where effects which are non perturbative 



in h may be important and one should solve the theory exactly |14j . 

In this paper we adopt the 1/N approach to quantum gravity. We then consider the 
following quadratic gravitational Lagrangian in four dimensions 



+ C n = J d 4 x | - 2A + R + aR 2 + (3R pv R» v + WttGC™} , (1.1) 
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where we have made use of the Gauss-Bonnet invariant to eliminate the dependence on the 
Riemann tensor. 

The semiclassical field equations derived by extremizing the action I can be written as 
Rp, ~ \Rg»v + A^, + a MHp, + (3 (2) i^ = 87rG(T^) ren , (1.2) 

where 

is the renormalized expectation value of the quantized matter source and where 

WH^ = -2R^ + 2g„ v UR - ~ 9tlv R 2 + 2RR flp , (1.4) 

and 

{2) H^ U = -2R« . ua + URp, + -g^DR + 2R«R au - ^g^R a(3 R al3 . (1.5) 

The (renormalized) values of the coupling constants a and (3 come from the regularization 
process and depend on the number and types of fields present as well as on the method of 
regularization. However, since the higher derivative terms must have very little effect on the 
low-frequency domain (in order to avoid conflict with observations), is usual to consider them 
as free numerically small parameters of the order of l 2 { (or equivalently, h, for G — c — 1). 

It was observed in 0,0], and independently in |JJ, that one can develop a perturbative 
approach to that problem that allows one to write the perturbative field equations in terms 
of second order derivatives of the metric and of the energy-momentum tensor. The system 
is thus reduced to an effective Einsteinian one. This perturbative reduction procedure can 
be described as follows. In the perturbative approach (I| we consider coupling constants a 
and (3 to be of the same order of l 2 u and curvatures small small enough to ensure that the 
perturbative series makes sense 

|ai2| « 1 , \pR pv \ < 1 • (1.6) 

Note that these conditions arise naturally since the description provided by the semiclassical 
gravity breaks down in a domain where the curvature is large i. e. near \R^ urJ i\ ~ c 3 /(hG) = 
3.9 x 10 65 cm- 2 . 

Bearing these conditions in mind, we observe that, up to first order in a and (3, the two 
geometrical conserved tensors and ^'H^, can be computed in terms of any known 

classical zeroth order solution, gj^J, of Einstein field equations 

^-i^{5+^ = ?^7£), (1.7) 

where Tjjj) is the stress energy momentum tensor of some matter fields, which can generally 
include both the source of an ordinary matter field Ff£) (e.g. electromagnetic field, perfect 
fluid, very energetic point-like particle, etc), given by T° v (°) = T^gffj, F^)), and the source 
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of A" identical scalar fields not interacting with the precedent ordinary matter fields, given 
by (°) = T^(g<$, <jf® j ) with j = 1, N. All matter fields are then coupled to g$, so 
that their only relevant back-reaction effects are on the gravitational field. Note that since 
we are interested in issues related to the back-reaction problem, the dynamics of the matter 
fields will not concern us here. Our purpose is to solve the equations of motion for gravity 
and to find the back-reaction of the matter fields on the gravitational field. 

From the complete knowledge of the zeroth order classical system, i. e. (gffi, Ff£), (f)^ J ), 
one can write the 'reduced' semiclassical equations to first-order (in a and 0) as Einsteinian 
like ones 




with an effective source Tj£) cfT . 

Let us now consider the following convenient approximation where the renormalized 
energy-momentum tensor in the expression of T^) eff can be splitted into two additive parts, 

(^> rcn = T« + (AT* «) , (1.9) 

where the first term of the right hand side of this equation includes both T° v W and 
T^ u W, and the second term includes only the (state-dependent) expected value of 
the energy-momentum tensor of the A" identical free scalar fields, (f)^ \ (AT^ U W) = 

(AT^ig®,}® j )) = 0{l 2 pl ). A similar approximat ion can also be found in Ref. IHl, in 
the linearized version of the one-loop approximation to quantum gravity. However, we note 
that the approximation ( |1.9] ) is more accurately justified in the context of 1/N expansion, 
where at least we have the following advantages: 

• The scalar fields <j>^ J are free (i. e. non-self- interacting), so that only one-loop terms 
in these matter fields arise. The stress-energy momentum tensor of these fields is 
approximatated in such a way that it can conveniently be splitted into a classical part, 
T^ u W, and into a quantum one, (AT^ V W). 

• We can use the same argument as that used for the one-loop graviton to neglect the 
one-loop (order contribution of matter field operator Ff£), so that the ordinary 
matter, T° u ^> = T^gj^J, F^}), can always be only treated classically. 

• The state-independent (local) corrections are uniquely determined (up the ambiguity 
given by the value of the coupling constants a and (5) by the two terms ^H^ u = 
W H ]j,u{g^J) and ^'H^ = ^ H^(g^J), which can always be build up from the zeroth- 
order classical metric gffl. 

A difficulty still remains. The state-dependent part (AT^ U W) requires further and 
independent work to be determined and, although we consider here only very simple quantum 
matter fields, to compute this term it is in general a very difficult task. But things can 
extremely simplify when the quantum state is build up from the conformal vacuum and 
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fields and background metric are conformally invariant. In this situation, in fact, all the 
corrections can be expressed in term of local quantities. 

Thus, once Tj£) cS is determined, one can try to solve the resulting reduced Einstein 
equations for the metric gj^J, which in turn will contain the first order corrections in a and 
(3 (§]. Note that, formally, the iteration can be pursued up to the n-th order to find Tf$ e ^ 
and thus one might solve the n-th order Einstein type field equations to find gj™) [f|, but 
since for many purposes it is enough and since explicit calculations are very complicated at 
order higher than the first we will stop the process to this order. 

A useful procedure, even to the first order in the iteration, is to solve for only some of the 
components of the reduced field equations and to enforce the conservation of the effective 
energy-momentum tensor, i.e. 

Tf y W ;I/ = . (1.10) 

This is the simplest way to ensure the existence of the solutions of the form gj}J, without 
having to solve the full set of field equations. 

It is evident, from the structure of the process, that we can only generate perturbative 
developments of the solution, if it exists at all. The knowledge of g^J allow us to search for 
first quantum corrections to several physical quantities of interest. 

It is also evident that for vacuum classical solutions the contribution of Q-'HfPj and ^H^J 
terms vanishes and the only non-zero contribution will be given by the term (AT^ W). This 
happens, for example, in the spacetime of a Schwarzschild black hole. Note that from this, 
it is also clear now, that the first order solution to the exterior metric of a gauge string given 
in Ref. ||, Eq. (23), has the last line wrong. In fact it gives a 1/r 2 term outside the core 
where Tf$ = and which is not present in the general relativistic solution. It is part of the 
motivations of this paper to correct this expression. 

The rest of the paper is organized as follows: Sec. II A. deals with the constant density 
model for the core of straight local cosmic strings. We look, first of all, for static state- 
independent perturbative solutions to this problem and we find that they do not exist. 
We then study in Sec. II B. (local) non-perturbative corrections to the metric by a method 
developed in Ref. that makes use of the properties of the fourth order field equations under 
conformal transformations (more generally Legendre transformations). Sec. II C. analyses 
the inclusion of some particular (non-local) quantum state term in the energy-momentum 
tensor that renders possible the existence of perturbative corrections to the metric. 

In Sec. Ill we study the issue of the existence of equilibrium solutions with spherical 
symmetry. Here we also take a constant energy-density p while we allow the pressure p to 
be dependent on the radial coordinate r to satisfy the Oppenheimer-Volkov equation. We 
again show the non consistency of this model and argue that non-local corrections to the 
renormalized energy-momentum tensor hardly can change here the situation. We end the 
paper with a further discussion of the perturbative approach to this uniform-density models. 

II. COSMIC STRINGS WITH UNIFORM DENSITY CORES 

In Refs. |2]|| the above perturbative method have been mainly applied to cosmological 
models. The main shortcoming here is the fact that we cannot trust the results in the most 
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interesting regime; near the initial singularity. The situation is less dramatic when dealing 
with black holes @||, since one is mainly interested in quantities evaluated outside or at 
the event horizon. The singularity at r = being hidden to the outside world until the very 
last stages of the black hole evaporation. 

On the other hand, we expect the perturbative approach to be very precise in the whole 
spacetime generated by a cosmic string (see conditions ( |1.6|) ), since (at GUT scale and 
below) curvature is moderate inside the string and vanishes outside. 



A. Perturbative solution with local source 



Let us consider a straight (gauge) string lying along the z axis. We have then a cylindri- 
cally symmetric system, Lorentz invariant in the z direction. In this case, a general enough 
metric can written as 



ds 2 = e 2Hr)( dt 2 + dz2) + dr 2 + r 2 e Mr) d(j) 2 



(2.1) 



within the following Gaussian coordinates range: — oo < t < +00, —00 < z < +00, < r < 
r_, < < 2tt. 

The non-vanishing components of the Einstein tensor for this metric read 



G r r = b'(r) (2a'(r) + b\r)) 



(2.2) 



Gl = 3b'{rf + 2b"{r) 



(2.3) 



G\ = a\rf + a'(r) b'{r) + b\rf + a"(r) + b"(r) . 



(2.4) 



Given the energy-momentum tensor (effective or not) and equating it to the Einstein 
tensor (over 8n since we choose units such that G = c = 1) we get a set of equations 
whose solution determine the form of a(r) and b(r). These equations can be supplemented 
with the conservation law which in the case of an energy-momentum tensor diagonal and 
only dependent on the coordinate r (in its covariant-contravariant form) takes the following 
simplified form 



TI „ + 2b' (r) { TI - 7*) + a'(r) ( TI - T 



^ = 0. 



(2.5) 



The 0-0 component of Einstein equations can be reduced to a Riccati equation upon 
the substitution y{r)=b'{r). Thus, 



2y\r) + 3y(rf = 8tiT* 



(2.6) 



Once solved this equation for y(r) by making the sum G r r + Gt — 2G\ we also get a Riccati 



equation for a'(r)=w(r) 

w\r) + w(r) 2 = 8n 



Tl 



Ti+n 



y{r) 



(2.7) 
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These two last equations together with the conservation law (|2.5|) are an equivalent set 
to Einstein equations (|2.2| )- (|2.4j ) and have the advantage of being of lower differential order 
[bearing only first derivatives in variables y(r) and w(r).} 

For the sake of simplicity we will now study the effect of the state-independent (local) 
semiclassical corrections to the uniform density model for the interior structure of straight 
cosmic strings. Thus the only non-vanishing components of Tj£) are |^-[T0[| 



I? = T«" = -^je (r ,-r), (2.8) 

where tq is a constant that specifies the energy density, r s is the "radius of the string" , and 
stands for the step function. 

It is easy to verify that in this case Eqs. (j2.5|)-(f2~T|) are satisfied by IPHTT) 



6(r) = (2.9) 



sin (j^ ; r < r s 
exp[a(r)] = { '° (2.10) 

'i)cos(^) ; r >r tan(^) , 



where we have imposed the equality of the intrinsic metric and the extrinsic curvature |TT 
on the matching surface located at r_ = r s and r + = r$ tan(r s /ro). Thus r s /ro parametrize 
the family of exterior solutions while r gives the scale of the interior problem. 

Given the above expression for the metric and Eq. (|2.8| ) we can compute the local part 
of the effective energy-momentum tensor from Eq. (|1.8|). Doing so, we obtain for r < r s 



rpt eff _ rp Z cfr _ _ ^ i (2a + (3) , , 

1 ~ z ~ 8vrr n 2 + 8vrr n 4 ' 



1 *« = T *<B = J 2a +f) . (2.12) 
v 87rro 

Since the local part of effective energy-momentum tensor vanishes for r > r s , the general 
form of metric exterior to the string remain unchanged. Only the parameter r s and the 
constant part of are expected a priory to vary. We will find, however, a somewhat 
unexpected result when we study the internal structure. In fact, since cS = constant, Eq. 
~M can be easily integrated. The general solution reads 



2 

b(r) = — In 



cosh ( J 6irTt eff (r - CV 



+ C 2 . (2.13) 



The constants C\ and C% can be determined from the matching conditions at r s with 
the exterior metric. In particular one would obtain C\ = r s and C 2 = 0. Note that since we 
have reduced the original fourth order set of equations (|1.2|) to an effective one of second 
order (Eq. ( |1.8| )), the matching conditions will be that of ordinary General Relativity JTT]: 
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6(r) 



ext 



b(r) 



hit 



and air) 



ext 



a{r) 



int 



d r b(r) 



ext 



r+ 



d r b(r) 



hit 



and d r a{r) 



ext 



r+ 



d r a(r) 



hit 



(2.14) 



i. e. the metric and extrinsic curvature should match as the boundary is approached from 
each side. 

One has also to keep in mind that since we are considering a perturbative approach, only 

eff 



up to terms linear in e have to be retained in ( |2.13| ), i.e. 



b(r) 



{2a + g) 



r - df + C 2 



(2.15) 



On the other hand, the conservation law 
leads to first order to 



1 



b\r) = . 



applied to T; eS (see Eqs. (|2~TT|)-(|2~T^)) 



(2.16) 



These two last equations are obviously inconsistent unless 2a + (3 = 0, what, in our case, 
leads us back to Einstein theory. Equivalently, if b'(r) = 0, Eq. ( |2.3| ) is not satisfied for 
the non-vanishing source (|2.11|) . What fails? Well, let us recall what our hypothesis have 
been. We have supposed that the semiclassical corrections to the gravitational theory are 
perturbative and can be absorbed in an effective energy-momentum tensor (see Eq. (1.8)). 
We then supposed that the resulting metric could be described by Eq. (|2.1|), cylindrically 
symmetric, Lorentz invariant and static. On this last point let us note that observe in Eq. 
( [2.13| ), although we have equal and constant energy density and pressure in the ^-direction, 
as in the general relativistic model that tend to balance each other, it also appears now a 
constant, but different from zero, radial and angular component of the pressure that would 
destroy the string since they are unbalanced. 

But let us keep the condition of staticity and explore a possible alternative which could 
make stable this simple model. 



B. Non perturbative solution with local source 

A second possibility is that the effects of the higher order terms in the Lagrangian ( |1 . 1| ) 
lead to local non-perturbative corrections. In that case, the method stated in Sec. I would 
give no result since by its procedure can only get perturbative corrections. In Ref. M it was 
developed a method to reduce the order of theory from fourth to second by introducing two 
auxiliary fields: A scalar field x an d a massive spin two field ipfiv These fields added to the 
usual graviton and account for the extra degrees of freedom of higher order theories. This 
method, although in a first approximation allows to compute non-perturbative corrections 
to the general relativistic metric. In fact, in Ref. it was found that (for the sake of 
simplicity we take (3 = 0) the metric solution to the higher order problem is conformaly 
related to the general relativistic metric g^ u 

9% = a + x(r))gl, (2.17) 
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where the field \ satisfies the equation 



(□ -ml)x 



8ttG 



rp(Matter) 



m 



2 • 



Qa 



(2.18) 



with the □ operator taken with respect to g E 



fJtU 



The metric outside the core of the string can be easily found from the above formulae |7J 

X (r) = CK (m r) , (2.19) 

where mo is real from the no-tachyons constraints 

3a + p>0, /3<0, (2.20) 

and C is a constant to be determined by the matching conditions with the internal metric. 

To find the internal metric is a little more involved. We find the solution x( r ) = Xp{ r ) + 
Xh( r ) where Xp( r ) — |( r o m o) * s the particular solution to fl2.18|) , and 



Xh(r)" + — cot (—) Xh{r)' - m 2 x h {r) = . 
r Vr / 



(2.21) 



The solution to this equation can be written in terms of Hypergeometric functions JT2 



Xh(r) = CiF [a,b, l,cos 



+ V C os 2fc — 
— V2r 



r 

{a)k{b) k 



r 

2^ 



In 



cos 



k=l 

-2tP(k + 1) + 2^(1) 
where coefficients 



+ C 2 ^F[a,b, l,cos 2 
ip(a + k) - ip(a) + ijj(b + k) - ifj(b) 



\2r J 



(2.22) 



1 + 



1 _ 2m 

1 3a 



2ro 
3a 



(2.23) 

shows the essentially non-perturbative character of the solution. Note that this non per- 
turbative solution is always convergent (and without poles) in the range of values of the 
coordinate r and, thus, well behaved at r = 0, what is in agreement with the hypothesis 
made in Ref. [0. 

The constants C, C\ and C 2 have to be determined with the matching conditions 



X(r) 



cxt 



X(r) 



int 



and d r x(r) 



drX{r) 



int 



[2.2A) 



but since both these external and internal non-perturbative solutions are already of com- 
plicate forms we do not write, here, explicitly their values. 

We finally observe, that although, in Ref. 0, it was emphasized the idea that a "self- 
consistent" method should be the only valid one to find "physical" solutions of the semiclas- 
sical field equations, in Ref. Q (where we developed a closely related approach) we never 
considered the perturbative method as an absolute physical prescription for "throwing out" 
all the non-perturbative solutions. Counterexamples to Ref. claim can, in fact, be found 
in recent papers ||T3|.|l~ 
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C. Perturbative solution with non local source 



Other possibility we have to discuss is to include in the effective energy-momentum 
tensor a state-dependent (non-local) piece. The local part of the total mean value is the 
one we considered in (|1.8|) . 

It is not difficult to compute (AT^ IU ) outside the string since the space is locally flat but 
topologically and, thus, globally different from Minkowski spacetime. In fact this non-local 
part will be the whole correction since the local part vanishes. One can perform the full 



quantum field theoretical computation [|15H T7|1 or make simpler symmetry and conservation 



arguments [T3j to arrive to 



Ah 

(AT^) = — diag(l,l,l,-3), (2.25) 



where the constant A depends on the number of conformal massless free fields we consider 



and we will keep its value generic [In Refs. [15-T7] it was found that for a single scalar field 
A = (14407r 2 )- 1 [cos- 4 (r s /r ) - 1].] 

Taking expression (|2.25| ) for the source of the semiclassical Einstein equations outside 



the core of the string and integrating expressions fl2.6|) and (|2.7f ), one gets for the corrected 
metric up to linear terms in h 

dg 2 = ( 1 _ A7r?lA \ ( ..2 . , 2\ , , 2 . (, , IMA\ 



with coordinate range: — oo < t < +oo, — oo < z < +oo, r + < r < +oo, < < 2n, 
and where now the boundary reads r + = r tan(r s /r ) + 27iAhcot(r s /r )/r . The spacetime 
possess an r-dependent deficit angle 

Ah 

D = 87Tfi-167T 2 —{l-Afi) (2.27) 

Now, in the interior of the cosmic string (r < r s ), we shall reverse the problem. Using 
similar arguments as in Ref. |18|, we will find the form of the (AT^) inside the string such 
that a static solution for the (classical) uniform-density source exists. Note that although 



in any static spacetime there exists a wide class of Hadamard states |T9| for which this 
expectation value makes physical sense, in our approximation ( |1.9| ), we chose the N scalar 
fields to be in the static vacuum state (i. e. (AT^) = (0|X^|0)), so that these fields do not 
contribute to the classical matter source (|2.8| ). 

Let us first note that when we include the non-local part of the energy-momentum 
tensor in the source of Einstein equations, and consider first order correction in a, f3, and 
h, a formal solution for the metric coefficients in terms of the source can be easily found. In 
fact, from ( |2.6|) and (|2.13|) we obtain 



b(r) = C 2 - (2a 4 + /3) (r-C 1 ) 2 + 4 7 r f dr> f dr" {AT*) , (2.28) 



while from Eq. (|2.7|) we find 
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air 



In 





c 3 


sin ( — ] 


+ — - 


L Vr /_ 





C 4 r + 



(2a + /3)r 
2rl 



cot 



to f -% f dr " sm2 « AI ?> + <AT *> - 2<AT < 



^2.29) 



where here = stands for up to first order terms in a, (3 and h. And the constants Cj have 
to be determined upon matching this metric with the exterior one ( |2.26| ). 

Again we find that the conservation law fl2.5|) proves to be very useful in analyzing the 
different possibilities. Up to first order corrections and in the background of the constant 
density string (|2.11|) one gets 



(Ar;>, r + 



cot 



r 



(at;)-(atj) + 



b'{r) 



. 



(2.30) 



Two cases can be straightforwardly analyzed: 



a. Non-constant radial pressure: (AT r r ) jr ^ and (ATJT) = (ATf). [Note that this last 
equality of components of the energy-momentum tensor is already fulfilled (c.f. (|2.12|) )by 
the classical and local one-loop parts of the T^ u .] 

In this case we have G r r — G^ = 0. From Eqs. ( |2.2|) and ( |2.3|) we get an equation for b(r) 
which can be immediately integrated to give 



6(r) = -D t cos ( - ) + T> 2 

,r , 



[2.31] 



The constants D\ and D 2 can be determined by the matching conditions with the exterior 
metric ( |2.26| ) at r = r s . Doing so, we obtain 



D 2 



An Ah cot 3 (g) 



o sm 



2nAh 



cot- [ — 

^0 



2f-ot- 1 — 1-1 



(2.32) 



Now, from Eqs. Q2.6|) and (|2.12j) we find 

(2a + 13) 



(at;) = (at;) 



+ 



■irrr? 



i-KrK 



r 

cos | — 



(2.33) 



From Eq. (|2.29| ) we see that to find the metric coefficient a(r) we need also to know 
(AT*). We can find it making use of the expression for the trace anomaly of a massless, 
conformaly coupled field (see [p]] for a review on the subject) 

h 



(AT) 



2880tt 2 

h 

2880vr 2 



AC^xpC^ + B (r^W - i/? 2 ) + COR 



4+4 



hB 



14407r 2 r^ 



(2.34) 
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where A and B depend on the type and number of fields considered, and the last equality 
defined B. 

We then find the other two components of the non-local part of the energy-momentum 
tensor 



(Ar/) = (AT^) = -(AT)-(Ar;) . 



(2.35) 



We can now explicitly perform the integral in Eq. (|2.29 ) 
a(r) = In 



sin 



D 3 2 „ 

H h -L>\ COS 

r 3 



(2a + (3)r hBr 
DaTq — h 



2r 3 



7207rr|] 



r 

cot ( — 

7*0 



^2.36) 



where the constants D 3 and _D 4 can be found from the matching conditions with the exterior 
metric ( |2.26| ) at r = r s , 



2r D! 



-cos [ — 

r 



2 (2a + (3)r s hBr s 
D A r Q — h — — 



2r 2 



720vrr2 



r 

cot ( — 

r 



8irAh 



r 



■ cor 1 — 

To 



^2.37) 



£4 



+ 



(2a + /3)r s _ fr5r s 

2^ 7207rrf 
/ TiBr s (2a + f3)r 
V720vrr§ 



r 



2r* 



- sin" ( 
cot 



2T>i sin (^) i67rAn 



2r n 



cot^ [ - 
r 



^2.38) 



b. Constant radial pressure: Let us consider the complementary condition (AT*') ^ 
(AT?) in the simplifying case (AT r r ) ir = 0. 

From the conservation law ( |2.5| ) we find that to first order in h 

(2.39) 

Plugging this equation into fl2T6|) we obtain a first order differential equation for (AT*) 
which can be readily integrated to produce 



y(r)=b'(r) = -4vrr cot ( - ) ( ( AT;> - (ATJ) ) . 



, (AT r ) sin 



fe) (2a + /3) 



COS' 



(n>) 



COS' 



Co) 



6Tcr? 



r 

tan I — 

r 



^2.40) 



where Ei is a constant of integration. 

Now, with this expression we go back to Eq. ( |2.6D which allows us to find 



b(r) = E 2 -^±p^-Anrl{(AT^)\n 



cos 



r 



+ £1 In 



cos fe)+ sin fe) 



cos 



Gro) 



sm 



V2ro) 



(2.41) 
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where E\ and Ei can be completely determined by the matching conditions. 

From the expression for the trace anomaly ( [2. 34 ) we find the other two components of 
the non-local part of the energy-momentum tensor 



(AT/) = (AT*) = - [(AT) - (AT;) - {AT*) 



[2A2) 



We can now explicitly perform the integral in Eq. (|2.29|) 

(2a + P)r (2a + (3) 



a r 



In 



r 

sin | — 

r 



+ E 2 - 



2rl 



2rl 



In 



cos — 



E 3 
+ — 



hr 



E (2a + (3)r 

E 4 r — Sir 1 5 h 

2tq 3607rrQ 



r 

cot ( — 



-8nE 1 In 



cos (^)+ sm (2^) 



+ 



sin 



(2.43) 



The constants E 3 and £4 can be then found from the matching conditions with the exterior 
metric ( [2.26 ) at r = r s . 

To summarize, we have found the explicit form of the expectation value of the stress- 
energy-momentum tensor that allow a consistent static solution for the interior of the 
straight cosmic string. More explicitly, we have imposed the semiclassical field equations 
and solved them after applying the reduction of order procedure to find that the existence 
of a perturbative solution imply that the expectation value of the stress-energy-momentum 
tensor depends on a and (3. Now, we really do not expect the first order non-local part to 
depend on a and /3, since it should be computed on the zeroth order background. This can 
be explicitly verified by an independent computation of the renormalized energy-momentum 
tensor using field theory methods, but this is beyond the scope of the present work. The 
form of the expectation value, here, have been chosen ad hoc to satisfy the perturbative 
static field equations. In fact, the dependence on the coupling constants comes from the 
imposition of the matching conditions that allow to determine the integration constants left 
in the resolution of the reduced field equations ( 2.2S )-( 2.30| ). They contain terms that com- 
pensate for the local contributions to the effective energy-momentum tensor. Clearly, this 
is an inconsistency that leaves us with the only option of the non-perturbative dependence 
on a and /3 of the straight cosmic string solution. 



III. SPHERICAL STARS WITH CONSTANT ENERGY DENSITY 

In order to further understand this non existence of perturbative solutions and see how 
general this phenomenon might be, we will study another system with uniform-density, but 
this time possessing spherical symmetry instead of a cylindrical one. 

We will represent the spherically symmetric metric by 

The non-vanishing components of the Einstein tensor for the above metric are 
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r(l 



-2A\ 



(3.2) 



-2A\ 



(3-3) 



-2A 



l\2 



$'A' 



A/ 
r 



(3.4) 



^ — W • 



(3.5) 



The Einstein equation derived from equating ( |3.2| ) to T/ can be formally integrated 

8ttG 



-2A(r) 



1 + 



r' 2 Tt(r')dr' 



(3.6) 



where we have imposed the condition A(0) = anticipating that we are going to consider 
systems that are will be regular at the origin of coordinates. 

The difference G r r — G\ = 8ttG(T^ — T/) leads to a formal integral of $(r) 



$(r) = $(0) - A(r) + AnG J\'e 2A{r,) (T r r (r') - T*(r')) dr' . 



(3.7) 



Instead of using the third independent Einstein equation above, we will consider the 
conservation of the energy-momentum tensor equation (which can be derived from the three 
independent Einstein equations anyway). The only relevant conservation equation will be 
the radial one 



T r r , T + 2 -(T;.-T e e ) + (T;.-Tt) $' = 0. 



(3.8) 



Combining this last equation with G r r = 8ttGT£ to eliminate $' we obtain a constraint 
on the energy-momentum components 



2A 



2r 



1-e 



-2A 



+ 87rGTV 



(t; - t\ 



. 



(3.9) 



This is the, so called, Oppenheimer-Volkov equation. 

Let us next consider a static perfect fluid described by the following energy-momentum 
tensor 



T$ = -p(r) , T r r = T e e = T* = p(r) , 



(3.10) 



where p is the energy density and p the pressure of an infinitesimal volume of matter. 

An exact solution can be found 0] when p = constant. In fact, Eq. ( |3.6| ) can be trivially 
integrated to give 



-2A(r) 



8vrGr 2 



(3.11) 
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It is also easy to integrate the resulting Oppenheimer-Volkov equation Q3.9|) 

p(r) = p 

3p + 3p c -(p + 3p c )Jl- 8 -^ 



(3.12) 



where p c = p(0). The radial coordinate of the star R is defined as p(R) = 0, and from the 
above expression can be found to be given by 



R 



3 


1 P + Pc 


8nGp 


p + 3p c 



(3.13) 



and thus the total mass M of the star will be given by M = A%GpR? /3. 
To complete the solution, <£> can be found from Eq. ( |3.7| ) 



$(r) = In 



8nGpR 2 



8nGpr 2 



(3.14) 



Now, the above solution will represent our zeroth order solution in the perturbative 
scheme of Sec. I. In building up the effective energy-momentum tensor we first observe that 
the interior solution we have just described is conformally flat, i.e. its Weyl tensor vanishes. 
In this case tensors and are no longer independent but related to each other 

by ^H^ u = 3 ( 2 'Hfj, v and a new geometrical, object can be defined in a four dimensional 
spacetime, which is conserved only in the conformally flat case 



2 1 1 

^Hfiv = —R^Rux + -^RR/iu + -zR\<tR Xo ? 9hv — -rR 2 gtiu 



(3.15) 



For the background metric ( ^.11| ) and ( |3.14j ) one obtains 0] the first order correction 



(3) TTV 



h8<irG) 2 p(p + 2p)=-Ap(r) 

3 7 



(3.16) 



^Ht = -\{8*GYp*= l -Ap 

3 7 



(3.17) 



where 7 stands for the coupling constant associated to ^'H^ in the Lagrangian, playing a 
similar role to a and (3. 

We first observe that the ^H^ v correction can be written in a perfect fluid form with 
the above first order modifications to the pressure Ap(r) and energy density Ap (which is 
still constant throughout the interior!). However, the bad news for this model are that these 
expressions for p + Ap and p(r) + Ap(r) do not satisfy the Oppenheimer-Volkov constraint 
( p.9|) , which for a perfect fluid takes the following simplified form 



dp 
dr 



ip+p} 



,2A 



2r 



1-e 



-2A 



+ 8nGpr' 



(3.18) 



We observe that again we find an incompatibility when we suppose that perturbative 
corrections to the uniform-density model exist. The key point now have been that the per- 
turbative method determines the form of the corrected energy-density and the r-dependence 
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of the corrected pressure. We thus no longer have the freedom to chosse p(r) such to fit 
the Oppenheimer-Volkov equation. One can actually checkQ that neither the terms coming 
from (a + (3 / 2>)^ H ^ corrections are able to satisfy Eq. (|3.9|). In passing we note that this 
terms do not generate corrections of the form of a perfect fluid as did those generated by 

Here the non existence of equilibrium solution seems even more grave than in the case 
of cosmic strings since we do not expect the non-local terms in the renormalized energy- 
momentum tensor to play any relevant role. In fact, the state-dependent part of the vacuum 
expectation value of the source vanishes if we suppose that, in addition the classical fluid, 
only massless conformally invariant free fields (such as photons and massless neutrinos) are 
present. In this case, in fact, since we are in a conformally flat spacetime there is no particle 
production and the Green's functions can be found by making a conformal transformation 
to Minkowski spacetime, finding the appropriate Green's functions there, and transforming 
back to the original spacetime. More complicated coupling can indeed be considered, but, 
again, it remains the question of its physical plausibility. 

IV. CONCLUSION 

Throughout this paper, we have followed the suggestions coming from the derivation of 
an effective action in the 1/N approximation to quantum gravity in the large N limit. We 
have, thus, assumed that there can exist interesting regimes of applicability of a semiclassical 
theory of gravity. We estimated that in some physical systems the quantum nature of matter 
is significant compared to the quantum nature of gravity which might remain negligible. In 
fact, although the semiclassical approximation breaks down at the regime of high (Planck 
scale I p i) curvature, such as in the final stages of an evaporating black hole or at very early 
times in the evolution of Universe, the effective theory should be valid in many interesting 
cases, i. e. when the curvature approach, but always remains (significantly) less than 
the Planck scale. This is the case of cosmic string and spherical stars where, although 
semiclassical corrections are always small, it may happen that these corrections do not allow 
the interior gravitational field to be perturbative in the coupling constants a and (3. In fact, 
we have seen that if we consider in the back-reaction problem only local contributions to 
the effective source of the semiclassical equations for uniform-density models of the core 
of cosmic strings and relativistic stars, we obtain that the interior metric acquires only 
non-perturbative corrections in the coupling constants a and /3, even if the source depends 
linearly on them. This shows that one cannot always simply neglect them or truncate the 
solutions to the first perturbative order || . Note the different situation here with respect to 
the case of the "instability of the Minkowski space". While to render Minkowski stable in 
Ref. |§ it was considered the truncation of the solutions to the perturbative ones and them to 
first order; here, we have that such static perturbations do not exist for the uniform-density 
models with cylindrical or spherical symmetry. 



1 These results have been verified by use of the program of analytic manipulation "Cartan", work- 
ing within "Mathematica" [21]. 
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We here recall that in the 1/N approach to the effective semiclassical action it is consis- 
tent to consider solutions beyond the first perturbative order p2|23|. In fact, in the case of 



large N limit, Eqs. ( |1.2| ) are exact, and therefore it make perfectly sense to look at solutions 
which are either perturbatively expandable in powers of a and (3 (to all the orders), as we 
did for the charged black hole case in Ref. || or not-perturbative in those parameters, as 
in the present paper. In this case, it makes also sense to look for exact solutions to the field 
equations (See Ref. |TB|) since they are obtained from the 1/N approximation to quantum 
gravity and they are exact to the leading order (in the case of large N limit). Note also that 
in the full 1/N approximation self-consistency only requires that these solutions should be 
of the same order of the 1/N leading order field equations. 

One can then consider the contribution of perturbative, but very particular non-local 
terms. This method (plus considerations of symmetry and the expression of the trace 
anomaly) allowed us to determine the state-dependent terms in the renormalized energy- 
momentum tensor that make the static perturbative solution to exist. What these states 
correspond to goes beyond the scope of this paper and should be determined independently 
by use of the standard methods of quantum field theory in curved spacetime. However, one 
can argue that this states appear to depend on a and f3 already at the first non-trivial order, 
while we would expect they only to depend on the zeroth order background metric. 

Finally, we know that the uniform-density model is a first approximation to the more 
involved situation inside the core of strings and neutron stars. In such case of energy density 
and pressure dependent on r the analysis should be possibly made numerically. 
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